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Abstract 

We compute the exponent 7 for self-avoiding walks in three dimensions. We 
get 7 = 1.1575 ± 0.0006 in agreement with renormalization-group predictions. 
Earlier Monte Carlo and exact-enumeration determinations are now seen to be 
biased by corrections to scaling. 



The self-avoiding walk (SAW) is a well-known lattice model of a polymer molecule 
in a good solvent; moreover, because of its semplicity, it is an important test-case in 
the theory of critical phenomena. Much work has been done in computing critical 
exponents by a variety of theoretical approaches - Monte Carlo (MC), exact enu- 
meration/extrapolation (EE), and renormalization group (RG) - with the aim of 
comparing these determinations with each other and with the experimental results. 
Small but persistent discrepancies have emerged between the predictions from dif- 
ferent theoretical approaches. For example, extensive studies have been done on the 
exponent v which controls the critical behaviour of the length scales. Early MC sim- 
ulations and EE studies |IJ in the '70s predicted v = 3/5 in agreement with the Flory 
theory. When the length of the walks which were simulated was increased [§, |3|, [| 
the value of v decreased to v — 0.592 ± 0.002. This value was also supported by 
extended EE studies which provided an identical estimate ||. On the other hand, 
field-theoretic RG computations || [1], ||, ^ [T(J persistently gave v = 0.588 ± 0.001 
or even slightly lower. The discrepancy was clarified when the MC studies were 
extended to much longer walks: because of strong corrections with non-analytic ex- 
ponent A fa 0.5, the asymptotic regime is reached only for very long chains, and the 
results from shorter chains are systematically biased upward [IT, [T2], [T3|, |TJJ. A simu- 



lation |14| with walks of length up to N = 80000, using a data analysis taking careful 
account of the corrections to scaling, gave v = 0.5877±0.0006 in good agreement with 
the renormalization-group estimates. Also universality is well satisfied: a simulation 
|I3| in a slightly different geometry provided v = 0.5867 ± 0.0013 (68% confidence 
limit) while a recent high-statistics simulation [T3] for the off-lattice Kratky-Porod 



model with excluded volume gave v = 0.5880 ± 0.0018. One may ask if the same 
phenomenon occurs for the other critical exponents. For the exponent 7, there are 
indeed significant discrepancies in the existing theoretical predictions. MC and EE 
studies provide the estimates 

1.161 ± 0.001 EE, Ref. [5] 

7 = { 1.1619 ±0.0001 EE, Ref. fT?| (1) 



1.1608 ± 0.0003 MC, Ref. p] 



On the other hand the e-expansion predicts a lower value J7|: 1.157 ^ 7 ^ 1.160. 
Indeed a Borel-type resummation gives 7 = 1.160 ± 0.004 and 7 = 1.157 ± 0.003 if 
one forces the expansion to reproduce the exact value in two dimensions 7 = ||. One 
can also use the scaling relation 7 = (2 — 77)1/ and the estimates for 77 and v: the 
unconstrained e-expansion gives v = 0.5885 ± 0.0025 and 77 = 0.031 ± 0.003 so that 
7 = 1.159 ± 0.005; using the exactly known values for d — 2, v — | and 77 = one 
gets v = 0.5880 ± 0.0015 and 77 = 0.0320 ± 0.0025 so that 7 = 1.1572 ± 0.0035. More 
controversial is the status of the expansions at fixed spatial dimension d = 3. In [[?) 7 
is estimated as 7 = 1.1615±0.0020 while in M the final estimate is 7 ~ 1.1613. These 
estimates depend crucially on the critical value of the renormalized coupling constant 
g*: in @, § the estimate g* = 1.421 ±0.008 is used , while [| uses g* = 1.422 ±0.008. 



However Nickel |1| has pointed out that the present estimates of g* could be slightly 
higher than the correct value due to a possible non-analyticity of the /3-function at 
~g* which is usually neglected in the standard analyses. A reanalysis of the series 
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[fTOfl indicates that g* could be as low as 1.39 and predicts 7 = 1.1569 + 0.10(g* — 
1.39) ± 0.0004. In this paper we present a high-precision MC study for the exponent 
7 using walks up to length N = 40000. Our results confirm the important role 
played by corrections to scaling; our final estimate 7 = 1.1575 ±0.0006 is significantly 
lower than previous MC and EE results but it is in perfect agreement with the RG 



predictions. Moreover it clearly supports the analysis by Nickel |18| , [10|] and the value 
~g* pa 1.39 for the renormalized coupling constant. In presence of strong corrections 
to scaling, in order to get a reliable estimate of the critical exponents, one needs to 
perform the simulation in the large- N regime. This is only possible if the algorithm 
at hand does not exhibit too strong a critical slowing-down. For the study of 7 
for SAWs on the lattice the best available algorithm is the join-and-cut algorithm 
||19|| : in two dimensions the autocorrelation time, expressed in CPU units, behaves 
approximately as jV~ L5 while in three dimensions the behaviour is expected to be 
iV~ L2 . The algorithm is thus nearly optimal. Another advantage of this algorithm is 
that it does not require the determination of the connectivity constant, at variance 
with more standard algorithms. The join-and-cut algorithm works in the unorthodox 
ensemble T^ tot consisting of all pairs of SAWs (each walk starts at the origin and ends 
anywhere) such that the total number of steps in the two walks is some fixed number 
N tot . Each pair in the ensemble is given equal weight: therefore the two walks are 
not interacting except for the constraint on the sum of their lengths. One sweep of 
the algorithm consists of two steps: 

1. Starting from a pair of walks (001,002), we update each of them independently 
using some ergodic fixed-length algorithm. We use the pivot algorithm [^(], [51], [3| 
which is the best available one for the ensemble of fixed-length walks with free 
endpoints. 

2. We perform a join-and-cut move: we concatenate the two walks U\ and 002 
forming a new (not necessarily self-avoiding) walk u conc ; then we cut uj conc at 
a random position creating two new walks 00' 1 and u' 2 - If u/i and u' 2 are self- 
avoiding we keep them; otherwise the move is rejected and we stay with u% and 

More details on the dynamical critical behaviour and on the implementation of this 
algorithm can be found in [HJ. Let us now discuss how the critical exponent 7 can be 



estimated from the Monte Carlo data produced by the join-and-cut algorithm. Let 
us start by noticing that the random variable Ni, the length of the first walk, has the 
distribution 

for 1 < Ni < Ntot — 1 where Z(N tot ) is the obvious normalization factor and is the 
number of walks going from the origin to any lattice point with N steps. For large 
iV we have 

c N Pa ao^N^ 1 (3) 

and thus the idea is to make inferences of 7 from the observed statistics of N\. Of 
course the problem is that (Q) is an asymptotic formula valid only in the large- N 
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regime. We will thus proceed in the following way: we will suppose that (H) is valid 
for all iV > N min for many increasing values of N min and correspondingly we will get 
estimates 7(iV m j n ); these quantities are effective exponents which depend on N min and 
which give correct estimates of 7 as N min and N tot go to infinity. The determination 
of 7 from the data is obtained using the maximum-likelihood method. We will present 
here only the results: for a detailed discussion we refer the reader to 
consider the function 

e {N )= (1 if N min < N x < N tot (4) 
mmV ; 10 otherwise V ; 

and let X be the random variable 



19]. Given N„ 



X = log[N 1 (N tot -N 1 )] . (5) 

Then define 

A {^minj - —Tq v- (P) 

\ N min I 

where the average ( ■ ) is taken in the ensemble T/v tot sampled by the join-and-cut 
algorithm. The quantity defined in @ is estimated in the usual way from the Monte 
Carlo data obtaining is this way X^c {N m i n ). Then j(N m i n ) is computed by solving 
the equation 

Xftr™(N min ) = [X\ t h^(N min ) (7) 
where, for every function of Ni, we define 

i f(N)] ( N ) - ggSr mw-'jNn - my-i 

The variance of ^(N min ) is then given by 

Var [W = Var ^Mc{N min )) 
{[X]X\ thA {N min )Y 

where [X;X] = [X 2 ] - [X] 2 . We must finally compute Var (Xff c s (N min )). As this 
quantity is defined as the ratio of two mean values (see formula (B|)) one must take 
into account the correlation between denominator and numerator. Here we have used 
the standard formula for the variance of a ratio (valid in the large-sample limit) 

Finally let us mention how to combine data from runs at different values of Ntot- 
The approach we use consists in analyzing the data separately for each N to t and then 
in constructing the usual weighted average of the resulting estimates ^(N min ) with 
weights inversely proportional to the estimated variances. Let us now discuss our 
results. We have performed high-statistics runs at N tot = 200, 2000, 20000 and 40000. 
The number of iterations is reported in Table |I|. The total simulation took 16 months 
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on an AlphaStation 600 Mod 5/266. In the same table we report also, for two different 
values of N min , the autocorrelation times for the observable 

V( AT \ - Nmin _ N m j n (~\~\\ 

1 mmj ~ (xe N -) ie N . ) 1 ' 

which according to (fl0|) controls the errors on 7. We use here a standard autocorrela- 
tion analysis J| with a self-consistent window of 15rj nt ,y, supplemented by an ad hoc 
prescription to take into account tails in the autocorrelation functions |14]]. The con- 



tribution of the tail to r int y amounts to approximately 20% for the two larger values of 
N tot , 2% for N tot = 2000, while for N tot = 200 we have been unable to identify any tail. 
The results are consistent with the expectation of r^.y ~ j\r~o.3 go keeping into 
account that the CPU-time per iteration increases as ~ jV~°- 9 ; we find that the CPU- 
time per independent walk increases roughly as iV~ L2 . Our raw data are reported in 
Table |2| and Table |3| where we give for various N min the estimates of X||gf (iV miri ). 
Let us now look at the results. Let us consider first N to t = 200. We see here that 
the estimates of 7 increase with N min and for N min = 50 they give 7 ~ 1.1605. This 



value is in agreement with Monte Carlo studies [16] performed in the same range of 
values of N and exact-enumeration studies [H, O. Consider now N tot = 2000. One 



can see that the estimates of 7 are lower and indicate 1.1580 < 7 < 1.1585: clearly 
the estimate at N tot = 200 was biased upward by the corrections to scaling. Let us 
now consider the results of Table ^| where we give the estimates of 7 coming from the 
weighted average of the results with N to t = 20000 and N to t = 40000. The estimates 
are extremely flat and agree within error bars from N m i n = 1 to N m i n = 8000; they 
indicate 1.1574 < 7 < 1.1578. Again the new estimate is lower than the previous 
ones. At least at N to t = 2000 there are still systematic deviations which are larger 
than the statistical error. Of course the same could be true for our results at the larger 
value of N tot . The most straightforward way to solve the question would be to have 
data at a higher value of N tot , say N tot = 4 ■ 10 5 , with comparable statistics. However 
this is not possible with current computer resources. We have thus simply tried to 
estimate the systematic bias by comparing the results with different N to t assuming 
that the systematic error vanishes as N tot (this is the behaviour one expects if the 
corrections to cn, formula (|3|), vanish as N~ l l 2 ). In this way we get an estimate of 
0.0002-0.0003 as the residual bias on the data at N tot = 40000. We have thus taken 
as our final estimate the value at N m i n = 2000: 

7 = 1.1575 ±0.0003 ±0.0003 (12) 

where the first error is the statistical one (68% confidence limits) while the second 
is a subjective estimate of the error due to the corrections to scaling. Let us finally 
discuss the determinations of the connectivity constant /i which is usually derived 
together with the estimate of 7. Restricting ourselves to the case of the cubic lattice 
- the discussion is completely analogous in the other cases — EE and MC studies 
give the following estimates for (3 C = l//x: 

0.213496 ± 0.000004 EE, Ref. [5] 

p e ={ 0.213499 ±0.000001 EE, Ref. (13) 



0.213492 ±0.000001 MC, Ref. 16 
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Ntot 


Niter 


T int,Y(l) 


T int, Y(1000) 


200 


5 • 10 8 


1.47970 ± 0.00055 




2000 


6.2- 10 8 


2.808 ±0.022 




20000 


10 8 


6.96 ±0.18 


9.35 ±0.21 


40000 


8.5- 10 8 


8.80 ±0.23 


10.78 ±0.24 



Table 1: Number of iterations and autocorrelation times for the various values of N tot . 



However these results have probably a large systematic error as the corresponding 
estimates of 7 are quite far from the correct value. To give an estimate of this effect 
we have considered the series 



S(S) 



uv=o 



A' 



(14) 



where cn is the number of walks of length N. If 5 = 7, S(5) has a simple pole as 
(3 = P c and thus, neglecting subleading corrections, (3 C can be determined by a simple 
Pade analysis. Using the values of reported in || [HJ, we find 

R = f0. 2134987 ±0.0000001 for 5 = 1.1619 , . 

^ c ~ I 0.2134811 ±0.0000031 for 5 = 1.1576 ^ ' 

where the "error" indicates the spread of the Pade approximants we have considered. 
Thus if use the value of 7 given by ]17|], see (0), we get a value for (3 C in agreement 



with their estimate (Tl3|). However if we use our estimate of 7, we get a (3 C which is 



off by five to twenty times the stated error bars. Of course we cannot trust (15) as 
a good estimate of (3 C as the analysis still does not take into account corrections to 
scaling; however this simple exercise shows that f3 c is lower than expected. It should 



also be noticed that this trend is already visible in (0) where (3 C obtained by MC 



is lower that the EE estimates. We thank Alan Sokal and Ettore Vicari for useful 
discussions. 
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N to t = 


200 




N tot = 2000 


■^V min 


ycens 


7 


N min 


vcens 


7 


l 


8.701168 ±0.000052 


1.15288 ±0.00011 


l 


13.298266 ±0.000068 


1.15782 ±0.00013 


10 


8.842084 ±0.000036 


1.15808 ±0.00021 


100 


13.452571 ± 0.000046 


1.15802 ±0.00028 


20 


8.943002 ± 0.000026 


1.15866 ±0.00034 


200 


13.551931 ±0.000033 


1.15811 ±0.00045 


30 


9.017434 ±0.000019 


1.15875 ± 0.00053 


300 


13.625478 ±0.000025 


1.15838 ±0.00071 


40 


9.074675 ±0.000014 


1.15999 ±0.00084 


400 


13.682081 ±0.000018 


1.1598 ±0.0011 


50 


9.119083 ±0.000010 


1.1605 ±0.0014 


500 


13.725970 ±0.000013 


1.1584 ±0.0019 



Table 2: Raw data and estimates of 7 for N tot = 200 and N tot = 2000. 



min 


vcens 
A MC 


(20000) 




vcens 
A MC 


(40000) 


7 


1 


17.90254 


± 





00026 


19.28850 


± 





00010 


1 


15763 


± 





00018 


200 


17.94330 


± 





00023 


19 


310211 


± 





000095 


1 


15762 


± 





00021 


400 


17.97715 


± 





00021 


19 


329406 


± 





000090 


1 


15758 


± 





00024 


600 


18.00687 


± 





00020 


19 


346998 


± 





000086 


1 


15765 


± 





00026 


800 


18.03376 


± 





00019 


19 


363292 


± 





000083 


1 


15760 


± 





00028 


1000 


18.05830 


± 





00017 


19 


378624 


± 





000080 


1 


15764 


± 





00030 


1200 


18.08095 


± 





00016 


19 


393077 


± 





000077 


1 


15754 


± 





00032 


1400 


18.10198 


± 





00015 


19 


406824 


± 





000074 


1 


15749 


± 





00034 


1600 


18.12172 


± 





00014 


19 


419918 


± 





000071 


1 


15744 


± 





00036 


1800 


18.14022 


± 





00013 


19 


432477 


± 





000069 


1 


15755 


± 





00038 


2000 


18.15760 


± 





00012 


19 


444488 


± 





000067 


1 


15749 


± 





00040 


2200 


18.17394 


± 





00012 


19 


456041 


± 





000065 


1 


15746 


± 





00042 


2400 


18.18930 


± 





00011 


19 


467169 


± 





000063 


1 


15739 


± 





00045 


2600 


18.20388 


± 





00011 


19 


477901 


± 





000061 


1 


15741 


± 





00047 


2800 


18.21771 


± 





00010 


19 


488255 


± 





000059 


1 


15739 


± 





00050 


3000 


18.230811 


± 





000093 


19 


498275 


± 





000058 


1 


15748 


± 





00052 


3500 


18.260909 


± 





000081 


19 


521898 


± 





000053 


1 


15753 


± 





00058 


4000 


18.287471 


± 





000069 


19 


543763 


± 





000049 


1 


15772 


± 





00066 


5000 


18.331323 


± 





000049 


19 


582986 


± 





000043 


1 


15749 


± 





00083 


6000 


18.364968 


± 





000034 


19 


617183 


± 





000037 


1 


1566 


± 





0010 


7000 


18.389939 


± 





000021 


19 


647298 


± 





000032 


1 


1586 


± 





0013 


8000 


18.407212 


± 





000011 


19 


673754 


± 





000027 


1 


1582 


± 





0017 



Table 3: Raw data for N tot = 20000 and N tot = 40000 and combined estimate of 7 
for various values of N min . 
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